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This  is  the  third  and  last  in  a series  of  reports  on 
high-speed  tensile -impact  work  performed  in  the  Pioneering 
Research  Laboratory  of  the  US  Army  Natick  Laboratories 
(now  the  US  Army  Natick  Research  and  Development  Command). 
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appeared  in  February  197^* 


The  first  report  dealt  with  the  basic  mathematical 
techniques  and  their  application  to  a linear  material. 

The  second  report  extended  the  method  to  non-linear 
materials  that  do  not  exhibit  creep  or  other  time  de- 
pendence . The  present  report  deals  with  non-linear 
materials  that  do  exhibit  time  dependence . Since  all  actual 
materials  exhibit  creep  and  relaxation,  it  is  in  the  present 
report  that  we  can  expect  to  get  the  best  test  of  the 
agreement  between  the  theory  and  the  actual  behavior  of 
materials . 


Body  armor  must  be  resistant  to  high-speed  impact  in 
order  to  stop  projectiles  and  shell  fragments.  Parachute 
straps  are  subjected  to  sudden  tensile  stresses  when  the 
parachute  opens  and  the  falling  load  is  suddenly  decelerated 
The  work  described  in  the  present  series  of  reports  re- 
presents progress  toward  an  understanding  of  tensile  impact 
but  leaves  much  still  to  be  done. 
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Fig.  1.  String  subjected  to  tensile  impact.  12 

Fig.  2.  Strain  pattern  for  nylon  yarn  calculated  by  15 

successive  substitution  (lines),  and  by  method- 
of-characteristics  ( points ) . 

Fig.  3.  Model  adopted  for  use  in  this  report.  18 

Fig.  4.  Lattice  of  points  at  which  calculations  are  26 

to  be  made . 

Fig.  5*  Fhoto  of  nylon  during  tensile  impact.  The  33 

data  obtained  from  this  photo  were  selected  for 


use  in  the  present  paper.  The  speed  of  impact 
was  90  m sec"-1-  and  the  successive  exposures  were 
made  at  the  rate  of  7500  sec”  . 


Fig.  6.  Various  stress-strain  curves  for  nylon;  A,  35 

semi-static;  B,  linear;  C and  D,  dynamic  but 
developed  with  different  assumptions. 

Fig.  7*  Comparison  of  results  of  Run  l4  (Curve  C,  38 

X = 1 ,14  = 60),  solid  line;  with  experimental 

results, /broken  line. 

Fig.  8.  Comparison  of  results  of  Run  l6  (Curve  D,  40 

X = l,f/  = 60),  solid  line;  with  experimental 
results/ broken  line. 

Fig.  9*  Comparison  of  results  of  Run  21  (Linear  4l 

curve,  R = 1.15,  X = 1,  LL  = 24), solid  line;  with 
experimental  results,  broken  line. 

Fig.  10.  Photograph  from  which  the  data  designated  NC-6  43 

were  obtained. 

Fig.  11.  Selected  data  from  the  experiment  NC-6.  Strain  44 

versus  time  after  impact.  Each  curve  refers  to 
a different  point  on  the  specimen. 
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Selected  data  from  the  experiment  NC-7*  Strain 
versus  time  after  impact.  Each  curve  refers  to  a 
different  point  on  the  specimen. 


A few  of  the  curves  calculated  in  our  attempts  to 
match  the  curves  in  Figs.  11  and  12. 


Fig.  14.  String  subjected  to  special  end  conditions 


Table  1.  Summary  of  Computer  Runs 


Glossary  of  Terns 


first  time-integral  of  the  repetitive  process, 
corresponding  to  the  lattice  point  in  the  ith 
column  and  the  row. 

second  time-integral  of  the  repetitive  process. 


= auxiliary  quantity  found  in  the  repetitive 
process . 

= velocity  of  propagation  of  strain  waves. 

= propagation  velocity  at  zero  strain. 


function  giving  the  initial  displacement  of 
each  particle  in  a string. 

function  giving  the  initial  velocity  of  each 
particle  in  a string. 

subscript  indicating  the  column  number  of  a 
lattice  point. 

subscript  indicating  the  row  number  of  a 
lattice  point. 

spring  constant  in  the  model  of  the  material. 

a repetition  index  (superscript)  in  the  re- 
petitive process . 

length  of  unstretched  sample . 

number  of  lattice -point  intervals  in  the 
length  L. 

2 

(c/c0)  , a given  function  of  the  strain  £ . 

displacement  of  a particle  of  material  from  its 
unstretched  position  x,  also  see  s’;  also,  in 
the  appendix,  a quantity  appearing  in  Laplace 
transforms . 

s/Lv0;  the  bar  is  dropped  after  the  quantity 
is  introduced. 

portion  of  the  displacement  s due  to  spring  1. 


portion  of  the  displacement  s due  to  spring  2 


displacement  of  the  particle  of  material  at  the 
lattice-point  ij  from  its  anstretched  position 

xi«5  • 


portion  of  sijdae  to  spring  2. 
time,  also  see  t. 

tc/L;  the  jar  is  dropped  after  the  quantity 
is  introduced. 


velocity  of  a particle  of  sample 


velocity  of  impacted  end  of  sample 


distance  of  a particle  of  material  from  the 
fixed  (left)  end  of  the  sample  in  the  an- 
stretched state,  also  see  x. 


x/L;  the  bar  is  dropped  after  the  quantity  is 
introduced . 

(appendix)  quantity  defined  by  eq  (10a);  also, 
a quantity  defined  by  eq  (37a). 

normal  strain  along  axis  of  string,  = 

portion  of  £ due  to  spring  1. 

portion  of  £■  due  to  spring  2. 


ratio  of  length  of  spring  2 to  length  of 
spring  1. 


damping  constant  of  the  dashpot  across  spring  2 


cnL,  the  normalized  value  of 


density  of  unstretched  material. 

normal  stress  across  spring  1 (total  stress 
along  axis  of  string). 

contribution  of  the  dashpot  to  . 

contribution  of  spring  2 to  (T*“  • 


RESPONSE  OF  POLYMERS  TO  TENSILE  IMPACT  3.  The 
Integral -Equation,  Successive -Substitution  Method 
Applied  to  Non-Linear  Materials  having  Time- 
Dependence  (Creep,  Relaxation) 


1 . Introduction 


The  present  report  is  the  third  of  a series  of  three 
that  describe  the  results  of  mathematical  and  experimental 
investigations  of  the  response  of  polymers  to  tensile 
impact.  Professor  Crout  is  responsible  for  the  mathe- 
matical developments,  Mr.  Pilsworth  is  responsible  for 
the  experimental  data,  and  all  three  authors  have  worked 
on  the  physics  of  the  problem  and  the  writing  of  the  reports. 
For  further  background  material  reference  1 should  be 
consulted;  it  is  the  first  of  the  present  series  and  is 
referred  to  as  Report  I.  Similarly  reference  2 is  the 
second  report  of  the  series  and  is  referred  to  as 
Report  II;  the  present  report  is  of  course  Report  III. 

The  work  was  undertaken  primarily  because  of  our  in- 
terest in  the  production  of  improved  body  armor.  Ex- 
perimental data  had  been  obtained  showing  the  propagation 
of  strain  pulses  in  nylon  string  after  tensile  impact  at 
moderately  high  speeds.  A theory  and  mathematical  tech- 
niques were  recognized  as  important  means  of  explaining  what 
was  observed  and  of  predicting  what  might  occur  under  other 
conditions . 


1.  P.  D.  Crout,  M.  N.  Pilsworth,  Jr.,  and  H.  J.  Hoge, 
"Response  of  Polymers  to  Tensile  Impact.  1.  An  integral- 
equation,  successive-substitution  solution  for  use  in  digital 
computers,  and  its  application  to  a linear  elastic  model", 
U.S.  Army  Natick  Laboratories,  Technical  Report  69-I8-FR 
(July  1968),  71  p;  AD-682714. 

2.  P.  D.  Crout,  M.  N.  Pilsworth,  Jr.,  and  H.  J.  Hoge, 
"Response  of  Polymeiato  Tensile  Impact.  2.  Extension  of 
the  Integral -equation,  successive-substitution  solution  to 
non-linear,  time -independent  materials",  U.S.  Army  Natick 
Laboratories,  Technical  Report  74-26-FR  (Feb  1974), 

76  p;  AD-7808H. 
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In  Report  I the  mathematical  method  was  developed  and  was 
applied  to  the  problem  of  strain  propagation  in  a linear 
material  having  no  creep  or  other  time  dependence.  Since 
the  solution  of  this  problem  is  known,  and  since  the  new 
method  gives  the  correct  solution,  its  validity  for  this 
case  has  been  established,  and  by  inference  the  new 
method  can  also  be  expected  to  give  correct  answers  to 
related  problems  for  which  satisfactory  solutions  have  not 
been  available . 


In  Report  II  the  iterative  method  of  Report  I was 
extended  to  non-linear  materials  and  was  applied  to  a set 
of  experimental  data  that  had  been  previously  analyzed  by 
the  method  of  characteristics.  Good  agreement  was  found 
between  the  present  (iterative)  method  and  the  method 
of  characteristics  at  low  and  intermediate  strains  but 
trouble  was  encountered  at  higher  strains  when  the  stress  - 
strain  curve  was  concave  upward. 


In  the  present  report  the  problem  of  a non-linear, 
time-dependent  material  is  treated  and  compared  with  ex- 
periment. Reports  I and  II  were  considered  as  setting  the 
stage  for  Report  III,  for  whereas  the  problems  treated  in 
I and  II  could  be  solved  by  other  methods,  there  is  no 
fully  satisfactory  method  of  computing  the  behavior  of 
non-linear,  time -dependent  materials.  The  results  being 
presented  in  the  present  report  do  not  clear  up  this 
situation  but  represent  some  progress  toward  solving  it. 
The  methods  now  being  presented  are  somewhat  limited  by 
the  tendency  toward  oscillatory  and  diverging  solutions. 
This  tendency  is  reduced,  as  we  expected,  in  the  presence 
of  time  dependency.  In  the  present  report,  the  greatest 
weakness  of  the  solutions  is  some  lack  of  uniqueness : 
there  is  too  much  latitude  in  the  choice  of  model  para- 
meters . A greater  range  of  experimental  data  would  un- 
doubtedly help  to  resolve  this  situation. 


. I 


The  experimental  data  with  which  we  compare  our 
results  has  mostly  been  obtained  on  specimens  50  cm  long, 
in  which  one  and  sometimes  two  reflections  of  the  strain 
pulse  occur  during  the  period  of  observation.  The  present 
report  contains  comparisons  of  computed  behavior  with  data 
in  which  reflections  occur,  and  it  also  contains  experimental 
data  taken  on  long  specimens  in  which  no  reflections 
occurred  during  the  period  of  observation.  The  no- 
reflection results  were  surprising  and  may  constitute  a 
major  contribution  of  the  present  report. 


. Summa: 


r of  Reports  I and  I 


-.f(5 O-t-tqW-tcrf J 


Equation  (2)  was  simplified  by  patting  f and  g=o  and 
introducing  three  new  normalized  variables : 


In  Report  I a new  method  of  handling  the  problem  of  strain- 
wave  propagation  was  developed,  in  which  the  relevant  partial 
differential  equations  are  transformed  into  integral  equations. 
The  integral  equations  are  then  solved  by  the  method  of 
successive  substitutions.  Referring  to  the  system  shown 
in  Fig.  1,  the  differential  equation 


in  which  v0  is  the  velocity  of  the  impacted  end  of  the 

string.  The  result  of  these  substitutions  (see  Report  I 
p 36)  was  to  reduce  eq  (2)  to  the  form 


was  derived  for  strain  propagation.  Here  x is  the  position 
of  a particle  of  material  in  the  unstrained  state,  s is 
the  displacement  of  the  same  particle  from  its  original 
position,  and  c is  the  wave  velocity.  For  ready  reference, 
these  and  other  definitions  are  included  in  a glossary  at 
the  beginning  of  this  report.  By  integrating  eq  (l) 
twice,  the  equation 


was  obtained;  this  equation  is  valid  when  c is  constant 
(a  linear  material  has  c constant).  The  function  f(x) 
specifies  the  initial  displacements  and  g(x)  specifies  the 
initial  velocities;  except  at  the  impacted  end  of  the 
string  these  functions  are  usually  zero  in  practical 
problems . 


=■  o. 


where  for  convenience  the  bars  have  been  dropped  from  the 
new  variables.  The  velocity  c does  not  appear  in  eq  (k); 
this  equation  refers,  in  fact,  to  the  problem  in  which 
L=l,  c=l,  and  a^=1.  When  a solution  of  eq  (U)  has  been 
obtained,  the  results  can  be  applied  to  actual  problems 
by  using  eq  (3).  The  method  devised  for  evaluating  the 
integral  in  eq  (h)  was  to  accept  a lattice  of  points  in 
the  (x,t)  plane,  to  specify  the  values  of  s along  three 
borders  of  this  lattice,  and  to  find  the  values  of  s at 
the  remaining  points  by  an  iterative  procedure.  To  find 
values  of  the  integrand  the  three-point  differentiation 
formula  was  used,  and  for  the  integration  the  trapezoidal 
formula  was  adopted.  Values  of  s were  calculated  one 
row  (one  value  of  t)  at  a time  by  iteration,  using  the 
border  values,  and  taking  all  other  starting  values  from 
the  previous  row.  When  the  specified  degree  of  con- 
vergence had  been  reached,  the  next  row  of  values  was 
calculated . 

Two  linear,  time -independent  problems  were  solved  in 
Report  I.  In  the  first  the  end  of  the  string  is  set  in 
motion  with  its  full  velocity  vc  at  t=o;  in  the  second 
the  velocity  rises  linearly  during  a finite  time -interval 
beginning  at  t=o.  The  solutions  found  by  iteration  were 
compared  with  the  exact  solutions,  which  are  known.  The 
agreement  was  good  except  where  the  exact  solution  had  a 
discontinuity  (abrupt  rise,  wave  front).  In  the  neighbor- 
hood cf  a discontinuity  the  values  (of  strain,  for 
example)  calculated  by  iteration  oscillated  above  and 
below  the  exact  values,  when  plotted  versus  x.  It  was 
clear  that  the  iterative  process  could  not  give  results 
that  reproduced  a discontinuity  exactly.  This  had  been 
anticipated.  In  adjusting  to  a discontinuity  the 
iterative  results  were  found  to  be  first  too  low,  then  too 
high,  and  so  on,  but  as  the  discontinuity  traveled  away 
from  any  particular  point,  the  amplitude  of  the  oscillations 

near  that  point  was  found  to  decrease . 

r 

The  second  of  the  problems  solved  in  Report  I,  in 
which  the  rise  in  velocity  occurred  during  a finite  time 
interval,  was  more  easily  handled  by  the  iterative  procedure 
than  the  problem  in  which  the  velocity  rise  was  in- 
stantaneous. This  was  of  course  expected.  There  was  still 
some  oscillation  of  the  iterative  values  above  and  below 
the  exact  solution  but  not  enough  to  destroy  the  usefulness 
of  the  solution.  The  conclusion  reached  in  Report  I was 
that  the  successive  substitution  procedure  was  satisfactory 
and  that  oscillations  relative  to  the  exact  solution  could 
be  troublesome  but  would  probably  not  be  serious . The 
greatest  chance  of  trouble  was  expected  to  be  associated 
with  calculations  based  on  a stress-strain  curve  that  was 
concave  upward;  in  such  cases  shock  waves  are  mathematically 
predicted,  for  example  by  the  method  of  characteristics . 
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In  Report  II  it  was  shown  that  a non-linear,  time- 
inde pendent  material  may  he  handled  in  a manner  similar  to 
a linear  material,  bat  with  eq  (4)  replaced  by 


(5) 


This  equation  differs2from  eq  ( h ) only  in  the-  presence 
of  the  factor  (c/c0)  • The  other  variables  in  the 

two  equations  are  the  same,  although  in  Report  II  and 
eq  (5)  the  bars  over  the  variables  have  not  been 
dropped.  The  propagation  velocity  c is  a function 
of  £ , and  c/c0  is  the  ratio  of  the  velocity  at  strain  £ 

to  the  velocity  at  strain  zero.  An  accepted  stress- 
strain  curve  permits  c/c0  to  be  found.  For  use  in  the 

repetitive  process  the  quantity  R = (c/cQ)  was  expressed 
in  a form  convenient  for  machine  calculation,  usually  as 
a polynomial  in  £ . The  function  R ( £ ) was  not  changed 
during  any  one  computational  run. 

Asa  test  of  the  usefulness  of  the  successive- 
substitution  method  for  non-linear,  time -independent 
systems,  a problem  was  selected  that  had  been  previously 
solved  by  the  method  of  characteristics . The  values 
of  R were  obtained  from  the  same  stress -strain  curve 
that  had  been  used  in  the  method-of-characteristics 
solution.  Good  agreement  was  obtained  between  the 
successive-substitution  solution  and  the  method-of- 
characteristics  solution,  from  the  time  of  impact  up 
until  the  time  when  the  method-of-characteristics  pre- 
dicted that  shock  waves  were  beginning  to  develop  and  the 
successive-substitution  method  began  to  give  serious 
oscillations.  The  two  sets  of  results  are  compared  in 
Fig.  2;  the  lines  show  strain  as  a function  of  position, 
calculated  by  successive  substitutions;  the  plotted 
points  are  corresponding  strains  calculated  by  the  method 
of  characteristics . All  the  points  and  the  corresponding 
curves  show  good  agreement  except  for  two  points  near  the 
top  of  the  figure.  These  points  and  the  corresponding 
curve  were  in  the  region  of  developing  shock  waves  and 
increasing  oscillations. 


STRAIN 


POSITION 


Strain  pattern  for  nylon  yarn  calculated  by 
successive  substitution  (lines),  and  by  method 
of -characteristics  (points). 


The  successive-substitution  solution  could  not  be  made 
to  converge  in  the  region  where  the  method  of  characteristics 
showed  the  development  of  shock  waves.  This  as  we  know 
is  the  region  where  the  stress -strain  curve  is  concave 
upward.  In  order  to  obtain  a solution  of  sorts,  the 
accepted  stress -strain  curve  was  modified  slightly  by 
replacing  the  final  part  of  the  concave -upward  section 
by  a straight  line;  this  avoided  divergence  of  the 
solution.  It  was  concluded  that  the  solution  does  not 
immediately  diverge  in  a region  where  the  stress -strain 
curve  is  concave  upward  but  diverges  when  the  computed  wave 
front  becomes  very  steep.  The  question  arose  as  to  whether 
the  shock  waves  indicated  by  the  method  of  characteristics 
and  the  oscillations  shown  by  the  successive-substitution 
method  were  simply  computational  results  not  found  in 
actual  materials.  This  was  a tempting  assumption,  since 
we  know  of  no  cases  where  a tensile  shock  wave  has  ever 
actually  been  observed.  However,  it  was  concluded  that 
the  possibility  of  actual  oscillations  in  strain  along 
the  length  of  a string  could  not  be  ruled  out . A con- 
siderable portion  of  Report  II  was  devoted  to  showing 
that  such  oscillations  can  in  fact  serve  as  a means  of 
conserving  energy.  In  this  connection  it  should  be 
remembered  that,  when  the  method  of  characteristics  is 
applied  in  a region  where  a shock  wave  develops,  the 
solution  satisfies  momentum  and  continuity  but  does  not 
conserve  energy.  In  using  it,  one  has  to  assume  an  un- 
specified mechanism  to  convert  the  excess  energy  into  heat. 

3 . Extension  of  the  Previous  Solution  to  Include  Time- 
Dependence 

The  results  obtained  previously  apply  to  either  a 
linear  material  (Report  I)  or  to  a non-linear  material 
(Report  II)  but  they  do  not  apply  to  a material  that  shows 
time -dependence  (creep  or  stress  relaxation).  Since  most 
materials  show  time -dependence  the  methods  of  Reports  I 
and  II  are  of  limited  applicability,  unless  experiments  are 
devised  that  take  place  so  rapidly  that  creep  and  re- 
laxation are  insignificant.  In  extending  the  previous 
solution  to  problems  involving  time -dependence,  we  wished 
to  choose  a model  for  our  material  that  could  adequately 
describe  the  time -dependence  and  at  the  same  time  introduce 
no  more  additional  mathematical  complications  than 
necessary. 


Choice  of  Model  System 


The  model  adopted  is  shown  in  Pig.  3«  It  consists  of 
a spring  and  dashpot  in  parallel,  with  a second  spring 
in  series  with  the  combination.  This  model  was  discussed 
briefly  in  Report  I.  It  has  been  widely  used.  Smith-^, 
for  example,  uses  it  in  a problem  involving  nylon  yarn. 
His  treatment  differs  from  ours,  however,  becau.se  his 
springs  are  linear  (stress  = a constant  x strain) 
whereas  our  springs  are  in  general  non-linear.  Both 
Smith's  dashpot  and  ours  are  linear;  that  is,  the  force 
supported  by  the  dashpot  is  proportional  to  the  rate  of 
change  of  strain. 


A three -element  model  was  accepted  as  probably  the 
simplest  that  could  show  the  desired  time  dependence. 
The  two-element  models  both  have  drawbacks . A spring 
and  dashpot  in  parallel  (Voigt  model)  obeys  the 
differential  equation . 


This  model  allows  creep  but  cannot  represent  stress  re 
taxation  at  constant  strain.  A spring  and  dashpot  in 
series  (Maxwell  model)  obeys  the  equation 


This  model  gives  an  incorrect  creep  behavior  because  the 
creep  continues  indefinitely  at  a constant  rate,  whereas 
the  materials  of  interest  to  us  show  a decreasing  creep 
rate.  The  accepted  model  (Fig.  3)  avoids  the  two  draw- 
backs just  mentioned,  but  of  course  its  validity  has  to 
be  established  by  experiment.  As  will  appear  later,  the 
model  has  been  found  adequate  in  the  sense  that  by 
adjusting  the  various  parameters  our  experimental  data 
could  usually  be  well  represented.  It  is  of  course 
possible  that  some  other  model  would  have  been  equally 
satisfactory.  It  is  even  possible  that  a simple  Voigt 
or  Maxwell  model  would  have  been  adequate  for  some  of 
our  data,  considering  that  we  were  dealing  with  high- 
speed impact  and  rather  short  intervals  of  time,  'but 
the  adequacy  of  these  simple  models  was  not  anticipated. 


3.  J.  C.  Smith,  "Wave  Propagation  in  a Three-Element 
Linear  Spring  and  Dashpot  Model  Filament",  J.  Appl . 
Fhys.  1697-704  (15  Mar  1966). 


The  analysis  of  the  model  shown  in  Fig.  3 was  carried 
oat  in  sach  a way  that  the  two  springs  could  obey  different 
stress -strain  carves,  bat  when  the  time  came  for  actual 
calculations  the  assumption  was  made  that  the  two  springs 
were  alike  except  in  "length".  The  length  of  spring  2 
was  taken  to  be  \ times  the  length  of  spring  1.  If 
the  same  stress  could  be  applied  to  both  springs,  the 
strains  would  obey  the  relation 


Of  course  the  stresses  across  the  two  springs  are  in 
general  unequal,  because  the  dashpot  across  spring  2 
supports  part  of  the  stress . The  ratio  of  spring 
lengths,  X > was  one  of  the  parameters  to  be  adjusted  in 
order  to  fit  the  experimental  data. 

As  a second  disposable  parameter  the  damping  constant  tx> 
of  the  dashpot  was  used.  In  the  simple  case  where  spring '2 
is  linear  with  spring-constant  K the  constant  u,  is 
simply  related  to  the  time  constant  of  the  system  con- 
sisting of  spring  2 and  the  dashpot  across  it;  we  have 


--Kr 


where  7*  is  the  time  constant. 


In  addition  to  assigning  the  values  of  \ and  JJL  we  also 
have  some  freedom  in  the  choice  of  the  stress -strain  curve 
used  for  both  springs . 


Mathematical  Analysis  of  the  Time -Dependent  Model. 

The  model  of  a time -dependent  material  shown  in  Fig.  3 
can  be  analyzed  by  making  certain  additions  and  changes 
in  the  methods  used  in  Reports  I and  II.  The  separate 
spring  1 is  equivalent  to  the  "springs"  handled  in 
Report  I<  The  spring  2,  with  the  dashpot  across  it,  forms 
a new  combination  that  has  not  been  previously  treated. 


More  specifically  we  shall  suppose  that  a cubic  centi- 
meter of  the  material  behaves  like  the  mechanical  system 
shown  in  Fig.  3.  Because  of  the  unit  length  and  cross- 
section  the  strains  may  be  considered  to  be  elongations, 
and  the  stresses  to  be  forces.  The  basic  equations  are 


a~  - 


+ c r. 
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(11) 

<r  « <r( <•,) 

(12) 

(rjsf 

(13) 

^ df 

(1U) 

We  thus  consider  the  elongation  £ to  be  composed  of  two 
parts,  one  of  which  depends  on  the  total  stress  alone, 
whereas  the  other  depends  also  on  the  time  (past  history). 
Similarly,  we  consider  the  force  q-  to  be  composed  of 
two  parts,  one  of  which  depends  on  alone,  whereas 
the  other  is  proportional  to  The  springs 

shown  in  Fig.  3 are  in  general  not  linear.  In  the  above 
equations  the  effect  of  inertia  is  not  included. 

In  the  case  of  our  string  the  above  relations  pertain 
to  each  point,  and  (14)  becomes 


<r,  * u iii.  • 

also,  we  regard  the  displacement  s to  be  composed  of  two 
component  displacements,  thus 


— A;  t 

(16) 

accordance  with  (ll) 

3 A*  « c 
<3%  * 

(17) 
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^ 
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Finally,  the  condition  for  dynamic  equilibrium  is  expressed 
by 


(See  Report  II,  eq  6). 


We  shall  consider  the  damping  coefficient  fJL  to  be  con 
stant,  which  corresponds  to  viscous  damping. 


Integrating  (15)  we  obtain 


The  problem  in  which  we  are  interested  is  one  in  which 
the  initial  distributions  of  displacement,  velocity,  and 
strain  are  zero.  In  this  case  (21)  and  (22)  become, 
respectively, 


rlf.VZ#* 


In  view  of  (12)  and  (l8),  (2k)  may  be  written 


c and  c0  hence  pertain  to  the  left  hand  spring  in  Fig.  3 
It  now  follows  that 


From  (27)  and  (28)  we  see  that 


(23)  now  "becomes 


or,  placing 


we  have 


■ I 


A possible  method  of  successive  substitutions  is  now 
given  by  the  following  equations,  which  are  used  in  the 
order  given,  and  in  which  the  bars  are  omitted,  it  being 
understood  that  the  various  quantities  are  normalized  as 
indicated  on  the  preceding  pages.  The  step  numbered) 

°r rak'™** ' are  indicated  by  superscripts.  Starting  with 
€ < ' and  we  obtain  ,(4+1)  and  c(*+0  by  the 

following  steps j ' c2 


We  note  that  in  this  successive  substitution  process  it 
is  not  required  that  a record  be  kept  of  s-^sg,  s or  € 

The  initial  and  boundary  conditions  are 


^(0,  ±>)  - o 
^ CX,  0)  - o 


0 < t , 


^ CX^O)-  o 0<:  % ^ ^ > 

ft****  oft . J <38) 

Since  these  conditions  are  on  s,  it  does  not  appear  that  we 
shall  be  able  to  confine  our  attention  to  the  strains 
alone.  Accordingly,  we  shall  place  (36)  in  the  form 

,Cfc+0  ft*  /*>  , . (39) 

0 *0 

Since  and  s2  are  both  positive,  it  follows  from  (38)  that 

^iCQf£)  ~ 0 OmJs  ^z(Q3±)  = 0 sUfUiMj.  0<£ 

A>,U,0)  - 0 OurU  ^i(X,O)=0 

From  the  second  of  these  equations  we  also  have 

6,(x,o)  = o e4(x,o)-oY  05X4  ' j 


also,  since 


I f 

'nOjt)  j± , 


and  since  <Tj  is  finite,  we  see  that  €aC/,ib)  cannot  be 
changed  abruptly  from  its  initial  value  of  zero  regardless 
of  what  motion  is  specified  for  the  right  end  of  the  string 
On  the  other  hand  €,C  hi)  may  change  abruptly.  Similarly, 
since 


we  see  that  for  any  fixed  value  of  x 
tinuous  function  of  t. 


We  shall  now  consider  a tentative  numerical  procedure 
for  applying  the  above  equations.  As  before,  we  shall  use 
a square  lattice,  as  indicated  in  Fig.  4.  The  various 
quantities  are  determined  one  row  at  a time;  hence,  we 
need  only  describe  how  to  obtain  data  for  the  row 
when  data  for  the  preceding  rows  are  available.  The 
various  steps  are  as  follows. 


l)  s is  known  exactly  at  i=o  and  i=n  for  all  values  of  j. 
As  a first  approximation  we  shall  take  the  values  of  s for 
the  ( j-1  row  as  those  for  the  j*'*1  row  when  0 ^ 

Similarly,  we  shall  take  the  values  of  sp  for  the  ( j-1  ) 

row  as  the  first  approximation  to  those  of  the  row. 
Subtracting  sp  ^ from  we  obtain  s, for  the  row 


2)  Knowing  S]_' for  the  row  we  determine  £f  for 

this  row;  and  hence,  using  (39)>  for  this  row. 

\ th  » CD 

3)  Knowing  (r  * for  the  j row,  is  determined  for 

this  row  using* (35)  • This  is  a subroutine  process  of 
successive  substitution. 


4)  Integrating 
obtained . 


numerically,  s, 


row  i3 


f0r  the  row  is  obtained  by  subtracting  Sp 
. We  note  that  s is  known  exactly  at  the  end 


from  s 
points • 


6)  We  now  proceed  as  with  step  2,  and  so  on 


. . 


r&n&rf 


In  the  various  numerical  differentiations  and  integrations 
appropriate  formulas  are  used.  To  fix  ideas  we  can 
tentatively  say  that  the  numerical  integrations  will  be 
carried  out  using  the  trapezoidal  formula  (ll6  Report  i); 
and  that  second  derivatives  will  be  obtained  using  the 
three  point  formula  (117  Report  i),  based  on  a 
quadratic  polynomial.  Except  at  end  points  first  de- 
rivatives will  be  obtained  using  the  three  point  formula 


based  on  a quadratic  polynomial;  at  the  end  points  they 
will  be  obtained  using  the  two  point  formula 


based  on  the  linear  polynomial,  the  upper  sign  being 
used  at  the  right  end,  and  the  lower  sign  being  used  a' 
the  left . 


Proceeding  along  the  lines  indicated  abcve,  the  set 
of  instructions  given  in  the  next  section  for  programming 
the  successive  repetition  process  for  use  in  connection 
with  a digital  computer  was  devised.  In  these  in- 
structions A is  a measure  of  the  first  time  integral  in 
(39)*  B is  a measure  of  the  second;  and  C is  a measure  of 
the  time  integral  in  (35)*  Instead  of  finishing  the 
calculation  of  for  each  round  of  the  main  successive 
repetition  process,  which  is  aimed  at  solving  (39) , she 
subroutine  indicated  in  3)  above  is  replaced  by  a single 
round  of  the  subroutine  process.  This  was  done  because 
it  was  felt  that  there  was  no  point  in  calculating  to 
an  accuracy  much  greater  than  during  the  course  of 
the  main  successive  repetition  process.  However,  if 
difficulty  is  encountered  with  convergence,  it  may  be 
necessary  to  include  the  full  subroutine  process.  Mean- 
while, until  such  difficulty  is  encountered  we  shall  use 
the  simpler  process  described  in  the  instructions . 


In  she  above  discussion  s' 


are  the  basic 


quantities  which  we  know  at  the  beginning  of  the  kch  round; 
and  s(k+l-)  and  s2(k+l)  are  the  quantities  which  we 

ultimately  calculate  during  the. course  of  the  round. 

If,  however,  we  use  s^'  andf^  ^'  as  the  basic  quantities, 

we  avoid  the  situation  wherein  £_  is  integrated  to  give 
S2;  whic*1  is  then  differentiated  to  give  £2  • result 

is  the  set  of  instructions  given  in  the  next  section. 

b , Description  of  the  Repetitive  Process  when  Applied 
to  the  Nonlinear  String  with  Creep  or  Relaxation 

Our  rroblein  Is  to  determine  the  vaLues  of  the 
quantities  s^  and  at  the  points  of  a square  lattice. 

i takes  the  values  0,  L,  2,  3>  >n;  and  j take3  the 

values  0,  l,  2,  3>  ’*  • n is  given;  hence  the  range  of 
values  taken  by  { is  bounded.  The  range  of  values  taken 
by  however,  has  no  upper  bound  except,  possibly,  one 
which  will  be  assigned  later  from  practical  con- 
siderations, Fig.  1. 

The  values  of  s±j  are  given  for  those  points  which  lie 
on  the  border,  thus 


r 0 


^yiO  ~ O 


= />2,3,  "'/H/  ] l(u6) 


>tyuj  j -0, 1,  2;  3.'". 


The  values  of  23ij  are  given  for  those  points  which  lie 
on  the  left  and  lower  parts  of  the  border,  thus 


2 A'lo'O')  i *0, J 


28 


The  valaes  of  s.^  and^sij  at  the  other  points  are  determined 

a row  at  a time that  is,  for  j fixed  and  i = 1,  2,  3» * * * > 

(n  - l)  for  sij,  and  n for  2sij  • The  valaes  of  j are  taken 

in  the  order  1,  2,  3>  *"•  In  view  of  this  it  will 
evidently  be  sufficient  to  describe  how  to  obtain  the  valaes 
of  sij  and  2sij  f°r  the  row  when  the  valaes  for  all 

of  the  preceding  rows  are  known. 


The  values  of  sj.j  and  2s ij  for  the  row  are  obtained 

as  the  limits  gotten,  respectively,  by  applying  repeatedly 
a process  shortly  to  be  described.  At  the  beginning 
of  the  Kth  repetition  we  have  the  quantities 


as  approximations  to  the  quantities 


respectively.  Here  k is  a superscript,  not  an  exponent. 
When  the  process  is  applied  starting  with  the  quantities 
(48),  the  quantities 


are  obtained  as  approximations  to  (49),  respectively. 

The  sequence  of  repetitions  is  started  by  taking  the 
values  of  s^ . and  2sij  for  the  (j-l)th  row  as  approx- 
J fh 

imations  to  the  values  for  the  j011  row,  respectively,  thus 


These  values  are,  of  course,  known.  The  sequence  of 
repetitions  is  terminated  when  for  sufficiently  large  k, 
the  output  (50)  duplicates  the  input  (48)  to  the  re- 
quired accuracy.  The  common  values  which  then  compose 
(48)  and  (50)  are  then  taken  as  the  final  result  (49). 
The  values  of  s^j  and  2s ij  which  compose  the  row  are 

thus  determined. 


Finally  we  shall  describe  the  process  by  which  we 
start  with  (48)  and  obtain  (50).  During  the  course  of 
the  calculations  we  compute  and  record,  in  addition  to 
the  quantities  si.  and  28iJ,  the  auxiliary  quantities 

5 2.^Aj  ) R*  j ) 2^ J;  A4J  , Bi-J  and  j • 

Along  the  lower  border,  where  j=o  these  quantities 
are  defined  to  be 


too  " 

0., 

£ 

r /o 

--0,  ,6lo--0, 

4 * 

) 

1 ~ Yik,0  ' 

: 0 

£ " 
z 00 

2^  10 

'0^2^  ;?.o 

**•> 

0 

R,o= 

/, 

l\c 

/ 

A*,©-  ^ 

/ 

'JoO 

— V) 

■ 0, 

= 

(") 

0 i 

<r 

?,0~ 

0 xo) 

, r;7>t 

’ 

-o’.; 

(0) 

0 , , 

°'»o=  0 

2.  (v  ) • 

- 0; 

0 - 

0 ^ (0  ) " 0 ^ * * 

* . . ,0 

•'  -0  • 
■ 0 

(0): 

0, 

A, 

0 • * 
> 

‘A,.., 

,0'  ' 

A 


0>2) 


B ,„  = 0,  IV- 0,  83o-0,‘-'‘,B„l.lyOji 

C0y°,  C2o-0,'”,  Cm-o'0*, 


30 


The  quantities  Rij>  Aj.j,  and  Bij  do  not  exist  at  the 

points  of  the  left  border,  where  i=o,  or  the  points  of 
the  right  border,  where  i=n.  As  is  the  case  with  sp j 

and  2s ij  the  auxiliary  quantities  t ^ > Rijj  ^2 

2&1J*  0^j>and  C>J  are  obtained  one  row  at  a 

time,  are  Rnown  in  the  first  (j-l)  rows,  and  are 
obtained  in  the  j^.row  as  the  limits  approached  by 
the  Quantities  ,6^5** , \ <^J  * a0"* -J  s 

3 ' y**  and  .(*D  J respectively,  during 
the  course^of  £he  successive  repetitions.  The  desired 
basic  process  is  now  described  by  the  relations 


5 • Application  of  the  Method  to  Experimental  Data  (with 
Reflections) 


In  order  to  test  the  model  selected  for  a time-de- 
pendent material  and  to  test  the  successive-substitution 
method  developed  in  the  preceding  mathematical  analyses, 
a set  of  experimental  data  obtained  on  a nylon  yarn  was 
selected.  This  set  of  data  was  originally  presented  and 
analyzed  by  Pilsworth  and  Hoge^ . It  was  used  in  Report  II 
as  the  basis  for  a comparison  between  the  results  of  a 
successive-substitution  analysis  and  the  best  of  the 
method-of-characteristics  analyses  given  ih  Ref.  4.  Mil 
of  the  data  were  obtained  from  the  high-speed  stroboscopic 
photograph  shown  in  Fig.  3 of  reference  4 , which  is 


4.  M.  N.  Pilsworth,  Jr.,  and  H.  J.  Hoge,  "Rate-De 
pendent  Response  of  Polymers  to  Tensile  Impact", 
Textile  Res.  J.  35,  129-39  (Feb  1965)- 


Fig.  5.  Hioto  of  nylon  daring  tensile  impact.  The  data  obtained  from  this  photo  were  selected  for  use  in  the 
present  paper.  The  speed  of  impact  was  9°  tn  sec-1  and  the  successive  exposures  were  made  at  the  rate 


reproduced  as  Fig.  5 of  the  present  report.  This  figure 
contains  6 usahLe  exposures  following  impact,  and  the 
strains  obtained  from  these  exposures  are  plotted  as 
the  broken  lines  of  Figs.  7,  8,  and  9>  where  they  are 
compared  with  results  from  three  different  computer 
runs . 

Each  attempt  to  calculate  the  behavior  of  the  nylon 
yarn  as  experimentally  observed  requires  a stress -strain 
curve  for  the  material  to  be  accepted;  it  also  requires 
values  to  be  chosen  for  the  parameters  A and  o,  . 

We  began  by  using  the  stress -strain  curve  that  mad 
previously  given  the  best  results . This  is  curve  C in 
Fig.  6;  it  is  also  curve  C in  Fig.  5 of  reference  4. 

This  curve  had  been  developed  by  successive  trials  in 
which  the  method  of  characteristics  was  used  and  it 
gave  a reasonably  good  fit  of  the  experimental  data, 
especially  the  early  data  obtained  before  the  second  or 
third  reflection  of  the  strain  pulse  had  occurred.  The 
values  of  A and  At  required  for  a calculation  were 
chosen  partly  by 'guesswork  and  partly  by  methods  similar 
to  those  described  in  the  appendix,  but  more  primitive. 

The  values  used  in  each  computer  run  are  given  in 
Table  1,  together  with  other  relevant  quantities 
describing  the  run.  After  each  run  had  been  plotted 
and  analyzed,  new  values  of  A > M.  or  other  relevant 
quantities  were  selected  for  use  'in  the  next  run. 

Since  no  table  of  computer  runs  was  included  either  in 
Report  I or  Report  II,  a complete  table  is  included  in 
the  present  report.  All  runs  that  yielded  useful  in- 
formation are  listed.  Omitted  runs  had  faulty  programs, 
diverged  very  quickly,  or  had  other  defects . 

Results  of  three  of  the  best  runs  are  shown  in 
the  present  report.  The  first  of  these,  Run  14,  is 
shown  in  Fig.  7,  where  strain  is  plotted  against  position, 
with  time  as  a parameter.  The  6 solid  lines  are  plotted 
from  6 rows  of  the  computer  print-out,  one  row  being 
chosen  to  correspond  to  each  of  the  6 usable  exposures 
in  Fig.  5*  The  common  time  of  each  exposure  and  the 
corresponding  row  of  the  print -out  is  indicated  near 
each  pair  of  curves  in  the  figure.  Good  agreement 
between  calculation  and  experiment  is  indicated  if  each 
pair  of  curves  nearly  coincides.  In  Fig.  7 the  agree- 
ment is  only  fair.  One  of  the  weaknesses  of  earlier 
calculations  has  been  taken  care  of.  Figure  2 shows 
regions  of  constant  strain  near  x = o for  two  curves 
( t = 0.479  and  O.613  )»  it  also  shows  such  regions 
near  x = 50  for  three  curves  ( t = 0.213,  0.346, 
and  0.479).  The  experimental  data  do  not  show  these 
flat  regions . Run  14  does  not  have  the  flat  regions 
and  in  this  respect  represents  an  advance  over  the 
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Table  1.  Summary  of  Computer  Runs 


All  significant  computer  runs  are  listed,  including 
those  made  for  the  present  report  and  also  those  made  for 
Reports  I and  II.  The  omitted  runs  had  fan1 ty  programs, 
diverged  at  an  early  stage,  or  were  considered  not  useful 
for  other  reasons . 


Run 

Stress  - 

Strain 

Relation 

R 

A 

h 

Notes 

1 

Lin 

_ 

_ 

_ 

a 

2 

ft 

- 

- 

- 

3-10 

C 
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- 

- 

b 

11 

ft 

- 

- 

- 

c 

12 

ft 

- 

1/3 

240 

d 

13 

tf 

_ 

1/3 

60 

14 

11 

- 

1 

60 

15 

Lin 

1 

1 

60 

g 

16 

D 

- 

1 . 

60 

6 

20 

Lin 

1 

1 

24 

21 

If 

1.15 

1 

24 

g 

22 

tf 

1.15 

2 

60 

23 

11 

O.85 

1 

60 

24 

tf 

0.85 

1 

24 

25 

ft 

0185 

2 

120 

26 

11 

O.85 

2 

60 

27 

L?„ 

O.85 

2 

24 

28 

11 

.70 

1 

60 

29 

11 

• 70 

1 

24 

30-33 

11 

1 

1 or  3 

7-5  or 
30 

e 

35 

11 

25/16 

1/9 

61.218 

36-41 

C 

— 

1 to 
25 

3-9  to 
30 

42 

43-49 

A 

Lin  (2) 

* 

5 

11.25 

f 

36 


Notes  for  Table  1 


a Runs  1-2,  treated  in  Report  1. 

b Runs  3"ii>  treated  in  Report  2. 

c Carve  C unmodified  (no  linear  final  section). 

d Runs  12-29,  treated  in  the  present  report;  they 
refer  to  short  specimens  with  strain  reflections . 

e Runs  30~33>  3^-^2,  treated  in  the  present  report, 
runs  without  reflections. 

f Runs  43-49,  treated  in  Appendix  2 of  Report  2; 
stress -strain  curve  consisted  of  2 linear  sections 
of  slopes  1:4. 

g Runs  15,  16,  21  gave  best  agreement  with  the  ex- 
perimental data  of  ail  runs  made,  with  21  probably 
the  best  of  the  3* 


Fig.  7*  Comparison  of  resalts  of  Run  14  (Carve  C, 

X = l,yi / = 60),  solid  line;  with  experimental 
results,  broken  line. 


calculations  shown  in  Pig.  2.  However,  the  slopes  of 
the  calculated  lines  in  Fig.  7 do  not  agree  with  the  data 
and  at  t = 0.613  the  slopes  even  have  the  wrong  sign. 
These  facts  made  it  seem  likely  that  the  propagation 
velocities  of  various  strain  increments  were  incorrect . 
Hence  we  concluded  that  the  faults  of  Run  14  could  not 
be  corrected  if  we  continued  to  use  curve  C as  our 
stress -strain  curve.  It  seemed  likely  that  curve  C had 
no  fundamental  significance,  that  it  was  not  a true  curve 
of  a time -independent  material  but  was  the  result  of 
trying  to  represent  time -dependent  data  with  a time- 
independent  model.  Presumably,  the  representation  was 
fairly  good  because  the  experimental  data  came  from 
a single  experiment  with  a single  time  scale.  If  a 
much  more  rapid  or  much  slower  experiment  were  to  be 
performed,  the  representation  based  on  curve  C would 
probably  be  worse • 


Assuming  that  curve  C was  not  of  primary  significance 
and  probably  described  the  behavior  of  the  material  only 
for  one  particular  speed  of  impact,  we  decided  to  try 
other  stress-strain  curves.  The  next  one  tried  was  a 
simple  linear  curve.  The  result  (Run  15)  was  judged  to 
be  better  than  any  previous  run.  In  Run  16,  the  results 
of  which  are  shown  in  Fig.  8,  a new  stress -strain  curve 
was  used.  This  curve  (curve  D of  Fig.  6)  was  developed 
by  the  same  methods  used  to  develop  curve  C;  these 
methods  are  described  in  reference  4 . Run  l6  was 

judged  to  represent  the  experimental  data  somewhat 
better  than  Run  15  and  hence  was  the  best  run  made  up 
to  that  date . 

After  3 unsuccessful  runs  (17,  18,  19),  a set  of 
10  runs  (20-29)  were  computed,  all  with  linear  stress - 
strain  curves . Various  values  of  A and  Ul  were  used, 
and  in  addition  the  value  of  R (=C^jc^as  assigned  values 
ranging  from  1.15  down  to  0.70.  For  a linear  curve,  £ 
is  constant,  of  course,  and  assigning  a value  of  I 

is  simply  a convenient  way  of  changing  the  value  of  c 
for  the  material.  Of  the  group  of  10  runs,  Run  21,  shown 
in  Fig.  9,  was  judged  to  represent  the  experimental  data 
best,  and  to  give  a slightly  better  fit  than  any  of  the 
earlier  runs. 

To  be  useful,  our  model  of  a time-dependent  material 
should  contain  only  a few  disposable  parameters,  and 
these  parameters  should  be  uniquely  determined  by  the 
experimental  data.  If  this  is  the  case  we  can  use  the 
model  with  some  confidence  to  predict  the  response  of 
the  material  under  other  experimental  conditions,  for 
example  when  the  impact  is  more  rapid  or  slower  than  the 
impacts  we  studied.  The  results  presented  in  Figs.  8 and 
9 show  that  our  model  is  not  unique;  we  were  able  to 


Position  cm 


Fig.  9.  Comparison  of  results  of  Run  21  (Linear 

curve,  R = 1.15,  A - i,yl»  - 24)  , solid  line;  with 
experimental  results,  hroken  line. 
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represent  the  experimental  data  almost  equally  well  with 
different  sets  of  parameters,  including  both  linear  and 
nonlinear  stress -strain  curves. 

One  way  of  improving  our  model  would  be  to  include 
more  than  one  time  constant  in  the  creep  and  relaxation 
mechanism.  Our  single  dashpot  has  essentially  only 
one  time  constant.  Som , authors  have  used  a wide 
spectrum  of  relaxation  times  to  explain  polymeric 
behavior.  However,  we  would  need  to  have  experimental 
data  covering  a wider  range  of  impact  velocities  before 
we  would  be  Justified  in  introducing  a larger  number  of 
elements  into  our  model. 

Since  observations  of  impact  at  substantially  higher 
velocities  would  have  required  major  changes  in  our 
experimental  equipment,  we  searched  for  other  ways  to 
cast  more  light  on  the  behavior  of  our  material.  The 
set  of  experimental  data  analyzed  above  were  taken  with 
a specimen  50  cm  long.  This  length  was  short  enough 
so  that  reflections  of  the  strain  pulse  occurred  during 
the  period  of  observation.  The  reflection  of  strains 
into  the  already-strained  material  made  the  data  more 
difficult  to  interpret . We  decided  to  perform  ex- 
periments in  which  reflections  did  not  occur.  Two 
such  experiments  are  described  in  the  following  section. 

6.  Experiments  Without  Reflections 

Specimens  of  nylon  yarn  400  cm  long  (8  times  the 
normal  length)  were  marked  in  the  usual  way  at  1-cm 
intervals  and  installed  in  the  usual  way.  The  camera 
was  aimed  at  the  end  of  the  specimen  where  impact 
occurred,  so  that  it  could  photograph  a length  of 
roughly  50  cm  of  the  yarn.  The  yarn  was  so  long  that 
the  photographing  was  completed  before  the  strain  pulse 
could  be  reflected  at  the  fixed  end  of  the  string  and 
be  propagated  back  into  the  field  of  the  camera.  After 
a considerable  amount  of  work  (synchronization  of  the 
impact  and  the  photographing  was  a problem),  two  good 
photos  of  impact  behavior  were  obtained.  The  data 
obtained  from  these  photos  have  been  designated  NC-6 
and  NC-7.  Fig.  10  shows  the  photo  from  which  NC-6  was 
obtained.  Figures  11  and  12  show  typical  curves 
plotted  from  each  set  of  data.  The  plotted  points  were 
derived  from  the  data  and  the  curves  were  faired  thru 
the  points  (they  are  not  computed  curves).  In  each 
figure  the  top  curve  refers  to  the  end  of  the  yarn 
attached  to  the  head  mass  and  the  bottom  curve  refers 
to  a point  near  the  other  edge  of  the  field  of  view  of 
the  camera;  the  middle  curve  refers  to  a point 
intermediate  between  the  other  two.  The  two  figures 
show  strain  as  a function  of  time  for  periods  of  about 


Fig.  11.  Selected  data  from  the  experiment  NC-6.  Strain 
versus  time  after  impact.  Bach  curve  refers  to 
a different  point  on  the  specimen. 


Fig.  12.  Selected  data  from  the  experiment  NC-7*  Strain 
versus  time  after  impact.  Each  curve  refers  to 
different  point  on  the  specimen. 


2 milliseconds  after  impact.  For  details  of  how  the  high- 
speed, stroboscopic  photographs  were  taken  and  how 
the  data  (positions  of  marks  on  the  strings)  were 
analyzed,  reference  1+  should  be  consulted. 

We  had  expected  all  curves  of  the  type  shown  in 
Figs.  11  and  12  to  show  an  initial  rapid  rise  corres- 
ponding to  the  arrival  of  the  strain  pulse,  followed  by 
a gradual  rise  corresponding  to  a further  stretching 
of  the  material.  Of  the  six  curves  plotted,  only  the 
top  curve  in  Fig.  12  shows  this  behavior.  The  top  curve 
in  Fig.  11  and  the  middle  curve  in  Fig.  12  show  essen- 
tially no  change  in  strain  after  the  passage  of  the 
initial  pulse.  The  three  remaining  curves  all  show  a 
decrease  of  strain  with  time.  Instead  of  continuing 
to  elongate  by  creep,  the  material  represented  by  these 
curves  appears  to  be  contracting  after  it  reaches 
maximum  strain.  There  are  wide  variations  in  strain 
at  different  points  on  the  specimen  at  a given  time  and 
also  at  a given  point  at  different  times.  For  example. 
Fig  11  at  t=2  msec  has  a strain  of  about  7$  near  the 
head  mass  and  a strain  of  less  than  3%  at  the  point 
33  cm  downstring.  Also,  at  the  33  cm  point  the  strain 
went  above  4$  at  0.6  msec  and  then  fell  to  below  3$ • 

The  strains  shown  in  Figs.  11  and  12  show  major 
departures  from  the  behavior  that  is  commonly  assumed, 
in  which  a strain  pulse  of  constant  height  travels 
along  the  specimen,  leaving  uniformly  strained  material 
behind  it.  The  possibility  that  the  data  are  in  error 
because  of  experimental  difficulties  has  been  con- 
sidered, but  we  are  not  aware  of  anything  that  could 
invalidate  these  two  experiments . The  fact  that  one  of 
the  curves  in  Figs . 11  and  12  slants  upward,  two  are 
horizontal  and  three  slant  downward  does  not  point 
toward  any  simple  type  of  experimental  error.  Wavelike 
irregularities  are  evident  in  some  of  the  curves,  es- 
pecially in  the  two  top  curves  that  describe  the  point 
nearest  the  head  mass . Such  waves  have  been  observed 
in  many  of  our  experiments  and  have  amplitudes  that 
sometimes  exceed  twice  the  scatter  of  the  data.  We 
believe  that  the  waves  actually  exist  but  can  only 
speculate  on  their  cause.  One  possible  cause  is  non- 
uniformity of  the  nylon  yarn. 


The  mathematically  simple  picture  of  impact  is  one 
in  which  the  head  mass  is  instantly  accelerated  to  some 
fixed  velocity  and  travels  with  constant  speed  there- 
after. Our  actual  experiments  differed  significantly 
from  this  simple  picture.  The  head  mass  had  a finite 
velocity-rise  time  of  roughly  0.3  msec,  as  can  be  seen 
from  Figs.  11  and  12  (the  head  mass  velocity  is  very 
nearly  the  same  as  the  particle  velocity  at  x=o). 


After  the  head  mass  reaches  its  maximum  speed  it 
gradually  slows  down  because  of  friction  and  because 
of  the  force  exerted  on  it  by  the  stretching  specimen, 
but  the  deceleration  is  too  small  to  be  very  noticeable 
in  Figs.  11  and  12.  The  acceleration  and  deceleration 
of  the  head  mass  complicate  the  situation  somewhat 
but  do  not  give  any  obvious  help  in  explaining  the 
curves  shown  in  Figs . 11  and  12 . 


The  data  taken  in  such  a way  as  to  avoid  strain- 
pulse  reflections  in  the  part  of  the  specimen  being 
observed  were  analyzed  by  the  same  methods  as  were  used 
in  section  5 • It  was  realized  that  if  our  model  of 
the  material  was  a good  one  it  should  predict  the  be- 
havior when  reflections  were  present  and  also  when  they 
were  not  present.  Thirteen  computer  runs  (30  - 42) 
were  made,  using  a variety  of  combinations  of  the 
constants  A and-  R and  various  types  of  stress  - 
strain  curves,  including  linear  "curves"  and  curve  C 
of  Fig.  6.  Since  the  specimen  length  was  now  8 times 
what  it  had  previously  been,  the  distance  between 
points  of  the  computation  lattice  was  increased  from 
1 cm  to  2 cm.  Even  so  the ^computer  time  per  run  was 
considerably  increased,  because  all  lattice  points  in 
each  row  had  to  be  computed,  although  the  strain  at 
most  of  them  had  not  begun  to  rise  above  zero. 


None  of  the  13  runs  gave  curves  that  looked  like 
those  in  Figs . 11  and  12 . A few  of  the  calculated 
curves  are  shown  in  Fig.  13*  In  only  one  run  did  we 
succeed  in  getting  the  strain  to  decrease  with  time; 
this  was  Run  32-  The  three  curves  for  different  values 
of  x differ  from  the  experimental  curves,  however,  in 
that  all  three  calculated  curves  approach  each  other 
and  become  nearly  coincident.  The  experimental  curves 
never  get  very  close  together  and  are  diverging  as  the 
experiment  ends . 


It  seems  best  not  to  try  to  propose  any  specific 
mechanism  for  the  behavior  shown  in  Figs.  11  and  12. 
The  results  were  unexpected  and  interesting;  they 
cast  some  doubt  on  the  validity  of  our  three-element 
model . However,  it  should  be  remembered  that  large 
amounts  of  energy  are  released  when  the  slider  is 
accelerated  by  the  rubber  bands.  Energy  could  travel 
more  rapidly  in  the  steel  tracks  and  possibly  in 
the  supporting  table  than  it  travels  in  the  nylon 
yarn.  It  seems  doubtful  that  any  mechanism  exists 
that  could  transfer  energy  to  the  fixed  end  of  the 
string  and  set  it  in  motion  before  the  pulse  in  the 
string  arrived,  but  some  such  mechanism  cannot  be 


A few  of  the  carves  calculated  in  our  attempts 
match  the  curves  in  Figs . 11  and  12 • 


completely  ruled  out . Further  experiments  are  needed, 
in  which  the  conditions  and  parameters  of  the  experiment 
are  varied. 

8.  Conclusions 


The  integral -equation,  successive-substitution 
method  of  analyzing  and  predicting  the  oehavior  of 
materials  under  tensile  impact  hac  been  extended  to 
include  time -dependency  (creep  and  relaxation).  A 
three-element  model,  consisting  of  a spring  and  dashpot 
in  parallel,  and  a second  spring  in  series  with  this 
combination  was  used.  The  behavior  of  short  (50  cm) 
specimens,  in  which  the  strain  pulse  experienced  re- 
flections, could  be  explained  by  this  model.  It  was 
in  fact  possible  to  describe  the  behavior  with  various 
models,  employing  different  stress -strain  curves 
and  other  different  model  parameters . Since  different 
sets  of  parameters  described  the  observed  behavior 
almost  equally  well,  no  unique  set  of  parameters  for 
the  material  could  be  chosen. 

When  long  specimens  were  used,  in  which  no  re- 
flections could  occur  during  the  time  of  the  ex- 
periment, no  satisfactory  model  of  the  process  was 
found.  The  material  was  observed  to  contract  after 
the  strain  pulse  had  passed  over  it  instead  of  re- 
maining at  constant  strain  or  continuing  to  elongate 
slightly.  A study  of  this  somewhat -unexpected  be- 
havior appears  to  be  the  most  fertile  area  for  future 
work,  when  and  if  impact  experiments  are  resumed. 
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The  purpose  of  this  appendix  is  to  record  certain 
ideas  that  were  being  considered  when  work  on  this 
project  was  discontinued.  These  have  to  do  with 
determining  ^ and  fJt  more  precisely.  At  present  we 
have  too  much  freedom  in  their  choice;  hence  we  shall 
consider  additional  restrictions  which  can  he  applied 
experimentally.  We  note  that  the  use  of  the  parameter^ 
implies  that  the  two  springs  are  of  the  same  nature, 
and  differ  only  in  length,  spring  2 being  \ times  as 
long  as  spring  1. 

Determination  of  A • Noting  Pig.  3 and  equations 
(LOT  - (lM,  we  see  that  under  slow  (quasistationary) 
loading  q - ~Q,  and 


gL<T~  _ rr-'/'r  \ At  I — <T"  (ty 

4,6  ' T 7a 


d€,  0~'(€, ) 


Under  fast  loading 


<f  = <r(6,)  4,  - [4  - i J. 


Initially  <T,=  0 ; hence  6=0,  and  (2a)  becomes 


(p  r;  (T~  (0)  £ 0 • (3a) 


It  follows  that 


77  - 


W&&AV  4j  — 0 « (4a) 
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Placing  the  initial  value  €/  = Oi.n  (la)  and  dividing 
into  (4a),  we  obtain  finally 


where ^ the  subscripts  0,  F,  and  S indicate  "initial", 
"fast",  and  "slow",  respectively.  This  gives  A after 
the  derivatives  on  the  right  hand  side  have  been 
determined  experimentally. 


Determination  of  . Let  us  suppose  that  the  springs 
are  linear,  and  have  spring  constants  K|_  and  Kg,  re- 
spectively. We  then  have 


e = +tz  , 


<r  = k,  e , = <r, + crj  = k2  e, 4- 


(6a) 


Starting  aj;  t=o  with  the  system  unstressed  lot  us 
maintain  t = / , then 


€,  t / 

K,  6,  ”^2^2.  * H*  ^2  “ 0» 


Taking  Laplace  transforms  we  have 

+ a.  Z. ' i - -x 


(8a) 


where  the  bars  indicate  transforms.  Solving  for 
obtain 


•/l  (Kz't’/fd,  + K,) 


where 


Noting  the  partial  fraction  expansions 


o(  A 0( 


and  using  a table  of  Laplace  transforms  we  find  that 
the  inverse  transform  of  ( 9a ) is 


Multiplying  by  to  obtain  CT , and  noting  (10a)  we 
obtain  finally 


Differentiating  with  respect  to  t gives 


■umj,  ^ 


from  which  we  see  that  C*  changes  exponentially  from 
one  value  to  another  with  a time  constant 


We  also  note  that  the  first  term  in  (14a)  is 
spring  constant  of  springs  1 and  2 in  series. 
is  determined  experimentally/  (15&)  gives  ^ 


Substituting  from  (L7&)  this  bacomes 


(21a)  may  be  written 


where 


and  —yi.  are  obtained  from  test  data. 

For  an  assumed  value  of  c - perhaps  the  initial  "wave 
velocity"—  yLt>/fic.a  and  i/\  can  be  determined  using 
least  squares  in  connection  with  (27a),  which  is 
written  for  several  points  in  the  xt  plane.  We 
note  that  (27a)  is  linear  in  f//  and  • 

An  approximation  to  c^  can  also  be  obtained  from 
(25a),  in  which  we  place 


wherein  K is  the  spring  constant  of  both  springs  in 
series  - obtained  by  a slow  (quasistationary)  test 
and  \ is  a tentative  value  which  is  later  improved 


1 

-2  Further  Constraint:: 

. U3e  of  Sinusoidal  Excitation. 

in  addition  to  tests  involving  transient  motion  it  is 
also  possible  to  devise  tests  involving  sinusoidal 
motion,  or  vibration.  For  example,  let  us  consider 
a string  of  length  L fastened  at  the  origin  and  ex- 
tending to  point  L on  the  x axis,  Fig.  lk.  At  the 
right  end  x-L  it  is  fastened  to  a device  capable  of 
exerting  tension,  both  constant  and  sinusoidal  of 
angular  frequency  (JJ  . First  we  prestretch  the  string 
so  that  it  has  a constant  strain  £ ; and  then  super- 

impose a small  sinusoidal  motion  on  this.  In  the 
treatment  which  follows  Kj  s where  6^  is 

the  initial  value  of  K4*K,Aand  5 ^ 

£ f 2 ^ Jfol  j and  4a.  are  the 'sinusoidal  parts Jof * 

the"7 s tress ea,  strains,  and  displacements  which  are 
the  excess  over  the  prestressed  values,  (incremental 
stresses,  strains,  and  displacements.) 

Using  complex  notation  to  take  care  of  the 
sinusoidally  varying  quantities,  as  is  commonly  done 
in  solving  alternating  current  network  and  transmission 
line  problems,  beam  vibration  problems,  etc.,  our  basic 
equations  become 


2 - d (F~ 

■f00  ^ ‘ cLx  ’ 

(29a) 

F - <r , +-  a\  5 

(30a) 

€ = &,  + , 

(31a) 

<r  = Ki 6,  } 

(32a) 

<rt  = fij  cx> 

(33a) 

= K.L  . 

(3^a) 

where  the  bars  indicate  the  complex  quantities  which 
correspond  to  the  various  sinusoidally  varying 
quantities,  and  j is  the  imaginary  unit.  It  follows 
that 

<T~  = : <T i "t 

= - K,  (35a) 

l = £,+  h 


Placing 


K/  (Kz  +j<v/4*)  ‘ 


(37a) 


(29a)  and  (36a)  give 


— 9-  = O 

eijc 


or 


£± 

ltc- 


■/*  (A  A* 


(38a) 


(39a) 


The  general  solution  of  this  differential 
equation  is 


Jb  - A Ajt/>bcOC  + BcobtkX.  (UOa) 
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K,t  Kj  +<JU3^ 


Substituting  (4ta.)  in  (i+9a)  we  obtain 


or,  noting  (37a) 


AiMjcf^  L 


At  the  end  x=L  this  becomes 


which  may  be  written 


where 


r 


1 


% 

Ife 

i? 


For  any  value  of  (JO  we  can  obtain (^experimentally . 

(53&)  then  gives  jQ  , after  which  (5ba.)  gives  o(  . 

"We  thus  know  ((  as  a function  of  (O  * This  knowledge 
is  used  in  connection  with  (37a)  in  order  to  obtain 
Ki  y Kg  y and  • 

In  so  doing  it  may  be  desirable  to  take  advantage 
of  the  fact  that  the  spring  contant  K for  both 
springs  in  series  is 


*/  & 

* ~ K,+  K*  » 


(55a) 


which  relation  can  be  used  to  reduce  the  number  of 
unknowns  by  one  if  we  are  willing  to  determine  K 
separately  by  obtaining  a slow  test  value  of  dL<r/*L6 
evaluated  at  £ a . 

In  any  case  we  can  write  (37a)  in  the  form 


-K,CKz+Ju>^)«*=  0,  (56a) 


insert  the  successive  pairs  of  values 
(tOtjCdi)  ^ t „ (CO^fOP^tained  by  experiment, 
(53a),  and  (5^-a);  and  aefermine  the  unknowns  so  that 


rtt  a 

X |pwt(MKa+Jc^O 


I z finjJttLsoi/m,  4 


(57a) 


* y(3  can  be  determined  from  (53a)  in  several  different 
manners,  among  which  are  the  following. 

1.  An  interation  procedure  can  probably  be  devised. 

2.  An  interpolation  procedure  can  certainly  be 
devised. 

3.  Using  a digital  computer  a table  of  values  of 
the  function  /3  tan  for  complex  arguments  can 
easily  be  constructed. 
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In  so  doing  a digital  computer  can  be  used  to  cover  the 
ranges  of  the  unknowns  - two  or  three  as  the  case  may  be, 
and  choose  the  best  combination  of  values . 


Even  without  a digital  computer  it  is  likely  that 
a least  square  procedure  can  be  devised. 

Another  possibility  is  to  write  (37a)  in  the  form 


K is  obtained  separately,  as  indicated  in  connection 
with  (55a);  then  since  and  X2  appear  linearly  in 
(59a),  a least  square  procedure  can  easily  be  devised 
to  determine  them  so  that 


Knowing  Xq,  X2,  and  K,  the  other  quantities  Kq,  K2, 
and  jjl0  follow  immediately . 

Since  K appears  linearly  in  (59&)  it  could,  if 
desired,  be  included  among  the  unknowns  instead  of 
being  determined  by  a separate  experiment. 


